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ABSTRACT 


S 


/ 

/ 

The  intensity  coherence  function  (ICF)  of  the  acoustic  wave  function  from  a  point 
source  is  derived  by  the  path-integral  technique  for  transmission  through  internal  waves 
in  the  presence  of  a  sound  channel.  Separations  in  time  are  emphasized,  although 
separations  in  transverse  horizontal  position,  vertical  position,  and  acoustic  frequency 
are  discussed.  Approximate  intensity  coherence  times,  lengths,  and  bandwidths  due  to 
internal-wave  fluctuations  are  derived.  Analytic  approximations  suitable  for  computer 
coding  are  presented  for  the  micropath  focussing  parameter  ty,  which  controls  the  devia¬ 
tion  of  higher  intensity  moments  from  the  Rayleigh-distribution  values. 


INTRODUCTION 


The  intensity  coherence  function  (ICF)  contains  important  statistical  information 
about  the  acoustic  field  that  has  traversed  a  medium  filled  with  random  fluctuations. 
We  show  in  this  paper  how  to  derive  a  number  of  results  about  the  ICF  which  were 
indicated  in  our  earlier  works  on  the  subject.1'3  In  particular,  we  show  how  to  derive 
the  ICF  from  the  path-integral  technique  for  transmission  in  the  presence  of  a  deter¬ 
ministic  sound  channel.  The  ICF  is  a  fourth  moment  of  the  field.  This  paper  follows 
logically  after  our  paper4  on  the  mutual  coherence  function  (MCF),  which  is  a  second 
moment.  We  also  present  some  analytic  approximations  to  the  micropath  focussing 
parameter  7,  which  controls  the  deviation  of  the  higher  moments  of  intensity  from  the 
Rayleigh-distribution  values.  These  approximations  are  suitable  for  computer  coding. 

We  derive  a  path-integral  expression  for  the  ICF  of  time,  and  then  give  explicit 
rules  for  calculating  it  in  the  special  case  of  internal-wave  medium  fluctuations.  We  are 
able  to  explicitly  evaluate  the  case  of  internal  waves  because  the  accepted  spectrum  of 
internal  waves  implies  a  nearly  quadratic  structure  function  and  solutions  are  known 
for  path  integrals  with  quadratic  actions.  In  our  results,  the  intensity  coherence  time, 
coherent  bandwidth,  and  coherence  lengths  are  determined  by  weighted  averages  of  the 
medium  fluctuation  spectra.  We  also  indicate  already-known  general  results  on  the  ICF 
of  space  and  frequency,  and  we  derive  a  new  result  for  the  relation  between  the  scales 
of  the  ICF  of  time  and  MCF  of  frequency  for  an  internal-wave  medium.  A  companion 
paper  presents  comparisons  of  these  theoretical  results  with  experiment.5 


I.  PATH-INTEGRAL  EXPRESSION  FOR  THE  ACOUSTIC  WAVEFUNC- 
TION 

We  begin  with  the  wave  equation  for  the  pressure  as  a  function  of  space  and  time 
in  the  presence  of  a  spatially  varying  wavespeed.  We  follow  the  notation  of  the  review 
article  by  Flatte':3  The  sound  speed  can  be  expressed  as 

C{  ?,  t)=C0\\+U0(s)+n{  x,  <)]  (1) 

where  C0  is  a  reference  sound  speed,  U0(z )  is  a  dimensionless  function  of  the  depth  z 
representing  the  deterministic  sound  channel,  and  n(  x,  t )  is  a  random,  zero-mean 


function  of  position  representing  the  effect  of  medium  fluctuations  such  as  internal 
waves.  The  wave  equation  for  an  acoustic  wave  is  unaffected  by  the  time  dependence  of 
fi  because  n  has  only  components  with  very  low  frequency. 

The  acoustic  pressure  ^  obeys  the  wave  equation 

dtl<f>-C2  W=0  (2) 

The  parabolic  approximation  consists  of  considering  solutions  in  which  waves  are 
travelling  only  at  small  angles  to  a  particular  direction;  in  the  ocean  this  direction  is  in 
the  horizontal,  labelled  by  x.  Thus  we  try 


4>=  exp[i(qx-<Tt)]il>(  ?  ,  0 


(3) 


where  q  and  a  are  the  wavenumber  and  frequency  of  an  acoustic  wave  travelling  along 
the  x  axis  at  speed  Ca\  that  is,  q^a/C^.  The  "reduced"  wavefunction  0  is  slowly 
varying  in  space  and  time  compared  with  q  and  a,  and  satisfies  a  parabolic  equation:6 

2iqdtrp=  {-dyy-d„+2q2[U0~Hi}  }  0  (4) 


Equation  (4)  is  a  Schrodinger  equation,  and  thus  its  solution  can  be  directly 
expressed  in  terms  of  a  Feynman  path  integral7 


0=/V  /  Dz 


exp{  iqS,,(  z  )-iq 


(5) 


where  the  path  integration  (indicated  by  D)  is  over  all  paths  z(x)  =  [y(x),  z(x)}  con¬ 
necting  the  source  to  the  receiver.  The  phase  associated  with  the  path  in  the  absence 
of  fluctuations  is 

R 

qS°=q  {  [  ^*y)2+  (6) 

and  N  is  a  normalization  factor  chosen  by  convention  so  that  0=1  for  (i=0. 


H.  PATH-INTEGRAL  EXPRESSIONS  FOR  THE  ICF 


The  ICF  is  <l(5)I(A)>  where  the  angle  brackets  indicate  averaging  over  the 
ensemble  of  random  n  functions.  The  ICF  measures  the  coherence  between  the  intensi¬ 
ties  at  two  different  points,  labelled  by  A  and  B.  These  points  may  be  separated  in 
space,  frequency,  time,  or  a  combination.  Consider  spatial  separations  that  lie  in  the 
( y,z )  plane.  The  ICF  is  a  fourth  moment  of  the  field;  we  designate  it  by: 


G4(B,A)  =  >  =  <|iV|4/ Dzy  Dz2  Dz3  DzA 


exp 


^S9(l)-iqS0(2)+iqS0(3)-tqSa(4)-iq  f  {/*(?  y)-fi(z2)-Hi(z3)- ft(z  <)}  dx 

0 


>  (7) 


where  paths  1  and  2  arrive  at  point  A  and  paths  3  and  4  at  point  B.  (All  paths  begin 
at  one  point:  the  source.)  The  ensemble  average  applies  only  to  the  n' s,  so  we  may  write 


G4(B,A)>=|N|4  /  DzyDz2Dz3D?i 


exp 


iqS0 (I)-  iqSt(2)+iqS,(Z)-  iqS0( 4)- 


where 


M=q2<  ( { 


> 


and  we  have  used  the  fact  that 

<exp  (ia)>  =  exp(—  -^-<a2>) 


(8) 


(9) 


(10) 


if  a  is  a  zero-mean,  Gaussian  random  variable,  such  as  any  combination  of  n’s  is 
assumed  to  be.  Even  if  a  is  not  a  Gaussian  random  variable,  (10)  can  still  be  true  if  the 
higher  moments  of  a  are  small. 

The  problem  of  finding  a  useful  result  for  the  ICF  now  reduces  to  evaluating  the 
quadruple  path  integral  in  (8).  We  first  expand  the  deterministic  sound  channel  function 
to  second  order  in  the  displacement  of  the  paths  away  from  the  equilibrium  ray.  That 
is,  we  define  zr(x)  as  the  function  that  satisfies  the  ray  equation  (in  the  parabolic 
approximation): 


4 


(ll) 

with  appropriate  boundary  conditions  determined  by  the  initial  and  final  points.  We 
expand  U0(zi)  around  the  point  2r(x).  Let  v,(x)  =  z,(x)  — zr  (x);  then 

&'.(*<)  »  +  -±-  U]  v?  (12) 

where  it  is  understood  that  U0  and  £/„"  are  evaluated  at  z,(x).  This  expression  will  be 
valid  as  long  as  the  effective  bundle  of  acoustic  energy  stays  well  confined  around  the 
unperturbed  ray,  cr(x).  If  more  than  one  solution  of  (11)  exists  (deterministic  mul¬ 
tipath),  then  we  treat  one  unperturbed  ray  at  a  time.  Addition  of  the  results  depends 
on  the  coherence  between  deterministic  rays,  which  can  often  be  well  estimated.  At 
caustics  our  method  breaks  down. 

The  quantities  v,(x)  can  be  thought  of  as  four  new  path  variables  defined  as  devia¬ 
tions  from  the  equilibrium  ray. 

We  can  make  a  transformation  of  the  four  path  variables  to  a  new  set  given  by8 

2  a  =  t>1  +  v2  +  y3  + 

i 

I 

2  §  =  yi  +y2  ~y3  -  yA 

2  7  =  v,  —  v2  -t>3  +  (13) 

- 

I 

l 

2  §  =  —  ^2  +  ^3  —  ^4 

The  a  path-variable  integration  can  be  carried  out  by  noting  that  VI  is  assumed 
to  be  independent  of  the  centroid  of  the  four  paths.  The  6  variable  is  then  forced  to 
be  zero  for  the  case  of  intensity  correlations. 

We  are  left  with  a  double  path  integral  over  §  and  7  for  the  intensity  correla¬ 
tion  at  two  nearby  points  a  distance  R  from  a  point  source.  The  two  points  A  and  B 
are  separated  transversely  by  zt  and  in  time  by  tt. 


5 


exp 


-iq 


(14) 


and 

R 

M=f[2d(  1,0)+2d(§,t)-d(§+1,t)-d(g-?,t)]dz  (15) 

where  d(  z,t)  is  the  phase  structure  function  density  defined  in  Esswein  and  Flatte'9 

The  dummy  path  variables  3  and  7  are  constrained  by  the  above  treatment  to  have 

boundary  conditions  given  by  j)  =  7  =0  at  x  =0  and  §  =  za,  7  =0  at  x=R .  The  0, 

and  7X  are  the  vertical  components  of  §  and  7  . 

At  this  point  it  is  worth  noting  that  a  number  of  endpoint  terms  that  resulted 
from  the  integrations  by  parts  have  been  subsumed  in  the  normalization,  which  is 
required  to  give  unity  at  zero  separation  of  the  two  receivers. 

m.  EVALUATION  OF  THE  MONOCHROMATIC  ICF  FOR  INTERNAL 
WAVES 

Under  the  Markov  approximation,  the  phase  structure  function  density  can  be 
expressed  as 

d{y ,  z,  0  =2q2  <^2>Lpf(y,  z,  t)  (16) 

where  <n“>  and  LP  have  been  defined  previously2  3  and  / (y ,  z,  t)  is  the  phase  corre¬ 
lation  function  (PCF)  defined  by  Esswein  and  Flatte'.9  The  PCF  has  been  evaluated  for 
internal  waves  by  Esswein  and  Flatte',9  using  a  combination  of  analytical  and  numerical 
techniques.  Since  y  and  z  are  both  functions  of  range,  x,  (16)  must  in  general  be 
evaluated  by  a  numerical  integration  code. 

At  small  separations,  approximations  to  the  PCF  are  possible.  For  example, 

/(0,0, (17) 
/( 0,  2,  0)  =-J-{*2}x2ln(x„/z) 


(181 


where  {<*^}  is  an  average  internal-wave  frequency  that  is  dependent  on  the  local  depth 
and  position  of  the  ray  zr(x),  and  likewise  the  quantities  {k2}  and  z0.  All  three  quanti¬ 
ties  are  evaluated  in  Esswein  and  Flatte'.9  Because  (18)  is  nearly  quadratic  the  loga¬ 
rithm  may  rather  accurately  be  replaced  by  a  constant.  One  of  our  problems  will  be  to 
establish  the  best  value  of  that  constant.  It  will  be  useful  to  remember  that  two  impor¬ 
tant  quantities  defining  the  type  of  fluctuations  are  <t>  and  A  where 

&  =q2  J  dx  <n2>LP  .  (19) 

and 

=q  J  dx  <H2>LP  {k$}  \g{x,  x)|  (20) 

_! 

where  g(x,x)  is  a  Green’s  function  defined  in  Flatte'.3  We  also  define  Lv  =  {it2}  2. 

Evaluations  of  <t>  and  A  for  some  particular  examples  are  given  by  Esswein  and 
Flatte'.10 

We  will  now  use  (16),  (17)  and  (18)  to  evaluate  GA.  We  need  to  note  that  the  hor¬ 
izontal  (y)  dependence  of  f(y,  z,  t)  is  very  weak;  therefore  we  will  approximate  the 
PCF  as  f(z,  t)  with  y  set  to  zero.  Then  (14)  becomes: 


G\  (2*  -  0=l'vl2  /  DdDiex p 


-iq  i[dz0d^-U';0-i]dx-^M 


(21) 


and 


M  =q2  J  <  n2>  Lp  [2/(7,  o)+2f(0,t)-f  (0+,,  t)-  f  (0  t) ]dx 


(22) 


with  endpoint  conditions  given  by  7(0)=7(/?)=0  and  0{O)  =0, 0(R)  =za  . 


IV.  EVALUATION  OF  THE  ICF  OF  TIME 

Consider  the  intensity  coherence  function  <I(f)I(0)>.  that  is,  set  za  =0  in  (21)  and 
(22).  It  is  not  difficult  to  show  that  the  main  contributions  to  the  double  path  integral 
come  from  two  regions,  one  where  0  <  Lv/<$  and  the  other  where  0  <  Ly/'P  The  two 

regions  are  independent,  so  first  consider  the  region  where  7  is  small.  In  that  case, 


7 


expansion  of  the  last  two  /’ s  in  M  yields 

R 

M  «  2 q2  /  <fi2>  LP  f(%0)dx  (23) 

o 

which  is  now  independent  of  t.  Since  we  normalize  to  the  answer  when  t  =0,  this  term 
gives  unity.  Now  the  term  where  0  is  small  yields 

R 

■V/«2  /  <n2>LP  [/(7, «)+/(£,  0-/(7,  t)]dx  (24) 

o 

which  is  separable  in  0  and  7.  To  evaluate  (21)  we  must  either  do  the  0  or  7  path 
integral  first.  It  is  of  interest  that  in  the  phase-screen  case,  one  does  the  7  integral 
first,  but  here  we  find  it  necessary  to  do  the  0  integral  first.8  That  is 


<I(t )I(0)  >=1  +N  J  D-y  exp 


<H*>LP 


•  j  D  0  exp 


-iq 


R  R 

J  [djd^-  U*0  l]dx- q2  f  <fi2> LPf  (0,  t)dx 
o' ■  O 


(25) 


We  know  from  the  phase-screen  example  that  the  0  integral  is  cut  off  by  the  /(/?,  t) 
term  in  the  exponential,  that  is  when 

R 

q2  J  <»2>  LPf(0,t)dx  «  1  (26) 

o 

but  this  is  approximately 

**(/-)2ln(Iv//J)-l  (27) 

LV 

Thus  the  important  contributions  come  when  0  «  Lv/$  as  expected,  and  the  logarithm 
may  be  replaced  by  In  <f>.  Now  the  0  integral  may  be  done  since  it  is  a  quadratic  path 
integral.  Let 

„  „  LP  ln$ 

H=q2<n2>-Z—  (28) 

Lv 

Completing  the  square  and  doing  the  0  integral,  which  is  now  subsumed  in  the  normali¬ 
zation,  we  are  left  with  the  7  integral: 


8 


i 


I  *  2 

<I(<)I(0)>=l+iV  /  D~i  expl  —  /  [<?„' 7  +  U,~1  ^dx 


^  ?2  /  </i2>77{^}  t2  ^  \n  )dx 

-  O  L,y  I 


(29) 


and  again  we  must  ask  what  value  of  qr  to  use  in  approximating  the  logarithm  by  a  con¬ 
stant.  Again  from  the  phase-screen  example  we  know  that  the  ~i  integral  is  cut  off  by 
the  geometric  term  (the  first  term)  in  the  exponent.  Appendix  A  shows  that  the 
appropriate  value  of  ~i  is  Lv  [2A4>(  ln<f>)'4)  in  order  to  make  the  geometric  term  of  order 
unity.  We  define 

1 


3,  = 


(30) 


2A<fr(ln  <*>)* 

and  we  note  that  we  are  going  to  replace  In  (L v/~i)  by  In  /?0. 

With  that  approximation,  the  qr  integral  is  a  quadratic  path  integral.  Although 
quadratic  path  integrals  can  be  done,  the  form  of  (28)  is  not  convenient  because  of  the 
appearance  of  the  function  H  as  a  coefficient  of  and  because  of  the  appearance  of 

Well-known  standard  techniques  can  be  employed  most  easily  if  we  return  to 

the  original  double  path  integral  (21),  with  our  new  understanding  that  M  can  be 
expressed  in  quadratic  form.  Remembering  that  we  are  setting  za  =0,  we  can  express 
the  ICF  in  terms  of  a  matrix  form  within  the  double  path  integral: 

<l(l)I(0)>  =  </2>  [l  +K  (<)]  (31) 


K  (t)=N  J  DdD~iex p 


■  1 1 !  (■>■« 
1  o 


0  1 
1  0 


2 i 

<7 


H  O 
O  J 


dx 


(3: 


where 


2>-JL{J>}  t2\n0, 

L  |/ 


(33) 
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3 


a 


where  N  is  determined  by  the  requirement  that  K(t)—*  1  as  t  — >0.  The  general  solu¬ 
tion  to  this  type  of  path  integral  is  given  by  solving  a  set  of  coupled  ordinary 
differential  equations  for  a  matrix  .V/: 


O  M= 0 


where  the  differential  matrix  operator  O  is  given  by  the  expression  inside  the  square 
brackets  in  (31),  the  initial  conditions  on  M  are  M(x  =  0)  =0,M'(x  =  0)=I,  and  the  solu¬ 


tion  is 


A'(0  = 


det  M0(R) 


det  M(R) 


where  A/„  is  the  solution  with  t  =  0. 


One  can  show  that  the  result  for  K(t)  is  independent  of  multiplying  the  H-J 
matrix  by  any  diagonal  matrix.  In  particular,  if  we  multiply  by 


c  =  —  I  us,  t 


we  find  our  ordinary  differential  equation  is 


„  [0  l)  „  LP  0  1 

(9,,+C.)  ,  o 


L/D  v  * 

— 7  Qf  ,  ■•>..  M  =  0 

1$  x  ) 


a  =ui  <(ln  <t>  ln/?„)2 


This  is  the  equation  that  is  used  in  evaluating  K(t). 

An  important  approximation  may  be  made  if  { w2}  is  not  a  function  of  x,  which 
would  happen  for  a  short,  horizontal  ray,  or  if  the  weighting  function  <p‘>Z.p/Z,v" 
were  strongly  peaked  around  upper  turning  points  so  that  only  { ur}  values  near 


horizontal  enter  in  a  significant  way.  In  that  case  the  two  equations  b*-r>.fnr  ri t ,  ,» 
to  each  other,  and  related  to  the  equation  for  the  micropath  bandwidth  fun ’  .■  n  O 
derived  in  Dashen  et  al.4  As  a  result,  one  can  calculate  K  from  the  equation 


K(t)=|Q(Aa)|2 


where  the  pseudofrequency  Act  is  given  by 


_A£ 1=^_  2{ 

a  ln4>  J* 


and  o  is  the  acoustic  frequency. 


V.  SMALL-TIME  APPROXIMATION  OF  THE  ICF 


Calculation  of  K(t)  is  cumbersome  and  lengthy.  It  is  useful  to  have  a  simpler 
approximation  based  on  the  limit  as  t  — *0.  It  can  be  shown  that  K(t)  approaches  unity 
in  a  quadratic  manner;  that  is,  there  exists  a  constant  i/'  such  that 


K(t)  -+\-va  t'2 


We  can  obtain  an  expression  for  i/'  immediately  in  the  case  that  (40)  is  valid,  (that 
is,  constant  {w2})  since  we  have  an  expression  for  Q(A<7)  as  Act— *0  .  From  (41)  and  Eq 
(65)  of  Dashen  et  al4  we  find 


i/'2  ,  of  o.  In 5, 

—  =4o*{^}  — 


If  { or}  is  not  constant  we  must  treat  (38)  when  a—*-0.  The  same  technique  that  led  tc 
Eq.  (67)  of  Dashen  et  al,4  only  in  matrix  form,  leads  to: 


tc'2  =72lnd(J!n4>  J  dz  <n2>  LP{  k{}  (  J2} 


/  dx'  <n2>LP{k-}  [s(x.x')]2 


where  it  is  understood  that  the  first  </t2>Ip(  kv] 2  is  evaluated  at  z  and  the  second 
at  x In  actual  numerical  evaluation  g(z.x')  can  be  broken  into  two  factors,  one  a 


function  of  x  and  the  other  a  function  of  x'  (with  different  functions  for  the  two  cases 


x>x');  this  can  speed  up  the  calculation  substantially. 


uv.vlv'. 


V  V  V 1 


11 


VI.  SEPARATIONS  IN  FREQUENCY  AND  SPACE 

We  collect  here  results  presented  in  other  papers  for  the  intensity  correlations 
between  receivers  separated  in  space  or  in  frequency: 


<iU)i(0)>= 


1  +exp 


<I2> 


(45) 


<I(A<7)I(0)>=[1 +|<3(Ao)|2]  <I2>  (46) 


where  D(z)  is  the  phase  structure  function  for  transverse  separation  z  .  The  use  of 

<I2>as  a  normalizing  factor  avoids  problems  with  the  first-order  correction  discussed 
in  Flatte',  Bernstein,  and  Dashen.11  The  above  results  imply  that  the  intensity  coher¬ 
ence  lengths  and  bandwidths  are  simply  ^2  smaller  than  the  acoustic-field  coherence 
scales  evaluated  in  Dashen  et  ai.4 


Vn.  HIGHER  ORDER  MOMENT  CORRECTIONS;  THE  MICRORAY 
FOCUSSING  PARAMETER  7 


We  present  here  some  analytical  approximations  we  have  used  to  calculate  -7, 
where  it  is  remembered  that  higher  intensity  moments  are  approximately  given  by11 


<!*>«  /V!<  />* 


2 


(47) 


In  terms  of  internal  waves  with  the  spectrum  assumed  in  previous  references,  the 
expression  for  7  is:11 

ft 
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where  kv  =nn(z)/neB,  k'v  =  irn(z')/n0B,  and  Mj  and  j*  are  constants  defined  in 
Flatte'.3 


In  the  case  of  no  sound  channel,  where  the  rays  are  straight  lines,  the  contribu¬ 
tions  to  the  x  integral  are  either  relatively  uniform  from  0  to  R  or  are  peaked  near  the 
end  points.  However,  for  curved  rays  the  contribution  may  peak  at  different  places 
along  the  ray.  To  avoid  doing  a  double-integral,  double-sum  problem,  we  need  approxi¬ 
mations 

First,  the  summation  over  J  is  similar  to  the  summation  necessary  to  calculate  the 
vertical  phase-structure-function  density,  and  yields 
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Then  we  may  express  P(j . jr)  much  more  simply  in  terms  of  two  functions  of  x. 
called  a  and  3: 

P(j  ,x)  =  (1  — cos  3j2)  exp  (-a;2)  (52) 


and 
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The  integral  expression  for  a  is  the  same  that  appears  as  the  inner  integral  in  r* 
and  1/^3 

We  still  have  two  integrals  and  a  sum.  The  sum  can  be  approximated  by  an 
integral  and  can  be  done  analytically,  if  a  is  small: 
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We  have  found  the  following  empirical  rules  for  evaluating  {P(j,x)}  :  First,  for 
very  small  a  <10-3  we  may  do  asymptotic  expansions  to  obtain: 

W.*)}=  7^-  (56) 

4  a 

Second  for  10-3  <a  <10-1  we  use  (55)  with  appropriate  expansions  of  the  Ei(x) 
function  for  small  argument.  We  retain  terms  to  x6. 

Third,  for  a>10-1  we  use  the  direct  sum  (50)  with  a  decreasing  number  of  terms  as 
a  becomes  large.  (The  cutoffs  we  have  used  for  1  through  6  terms  are  at  values  of  or 
equal  to  5.,  0.64,  0.33,  0.2,  0.13,  and  0.1.) 

Thus  we  have  reduced  our  problem  to  a  double  integral,  with  the  inner  integral 
necessary  for  the  calculations  of  several  interesting  quantities:  r0,  the  width  of  the 
micropath  bandwidth  function  Q;  u1,  the  time  scale  of  intensity  decorrelation;  and  now 
■y,  the  micropath  focussing  parameter  that  controls  the  higher  moments  of  intensity. 

Vm.  SUMMARY  AND  CONCLUSIONS 

The  derivation  of  the  intensity  coherence  function  (ICF)  of  time  by  the  path- 
integral  technique  has  been  given.  Allowance  for  a  deterministic  sound  channel  and  the 
presence  of  reasonable  inhomogeneity  and  anisotropy  in  the  fluctuation  field  has  been 
included.  The  ICF  has  been  evaluated  for  fluctuations  dominated  by  internal  waves, 
which  have  a  vertical  structure  function  that  is  nearly  quadratic.  Reasonably  simple 
expressions  in  terms  of  environmental  measurements  for  acoustic  intensity  coherence 
times  and  coherent  bandwidths  have  been  given. 

All  of  the  parameters  necessary  to  approximate  the  various  ICF’s  can  be 
evaluated  on  a  mini/micro  computer  in  a  few  minutes  for  a  typical  ocean-acoustic 
experiment. 

We  are  grateful  for  a  number  of  useful  conversations  with  Dennis  Creamer  and 
Frank  Henyey  of  CSND,  and  for  the  hospitality  of  the  La  Jolla  Institute,  where  some  of 
this  work  was  done.  This  work  was  supported  by  the  Office  of  Naval  Research  and  the 
Defense  Advanced  Research  Projects  Agency. 


Is  order  to  estimate  the  appropriate  value  of  7  to  insert  in  the  logarithm  in  (29), 
we  need  to  find  the  typical  order  of  magnitude  that  7  takes  on  within  the  path  integral. 
We  argue  that  this  typical  value  is  that  which  makes  the  first  integral  in  the  exponen¬ 
tial  of  order  unity.  (Larger  7  would  make  negligible  contribution  to  the  integral 
because  of  the  exponential  cutoff;  smaller  7  has  smaller  phase  space.  As  long  as  it  is 
small  enough,  the  second  term  in  the  exponential  is  less  important;  the  crossover  time 
where  the  second  term  will  dominate  occurs  at  A4*2  times  the  intensity  correlation 
time.) 

We  can  express  the  requirement  on  7  as 

~  q2J[0<f>ydx  «  1  (A.l) 

where 

i  _i 

0  =  H  2(5„+u;)H  2  (A.2) 

i 

d>  =H27  (A.3) 

Thus  O  is  an  operator  and  we  can  imagine  its  eigenvalues  and  eigenfunctions,  with 

tr  0-,=/  H  (x)g(x,x)  dx  =  £  —  «  —  (A.4) 

where  g{x,x)  is  the  Green’s  function  for  the  operator  dzx  +U^  defined  in  Flatte',3  and 
the  lowest  eigenvalue  is  C0. 

The  definitions3  of  the  diffraction  parameter  A  and  the  strength  parameter  can 
be  used  to  estimate  £e  knowing  H  from  (28): 

C^ssgA^ln#  (A. 5) 

Let  us  define  70  as  our  estimate  of  the  order  of  magnitude  of  7.  We  can  make  such  an 

estimate  in  many  different  ways;  we  choose  the  following: 


which  is  an  estimate  with  weighting  along  the  ray  determined  by  H(x).  Noting  that 


we  can  express  (A.l)  as 
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and  we  can  also  use  the  definition  of  the  strength  parameter  $  to  show  that 


It  follows  that 


~  2L„  A<£  (ln<I>)% 


which  provides  the  appropriate  expression  for  the  definition  of  0O  in  (30). 
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